Rotational optomechanical coupling of a spinning dielectric sphere 
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We formulate a non-relativistic Hamiltonian in order to describe how the rotational degrees of 
freedom of a dielectric sphere and quantized light fields are coupled. Such an interaction is shown 
to take a form of angular momentum coupling governed by the field angular momentum inside the 
dielectric. As a specific example, we show that the coupling due to a single whispering gallery mode 
can lead to precession dynamics and frequency shifts of light. 
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Cavity quantum optomechanics has been an active re- 
search area investigating quantum phenomena and ap- 
plications through the interaction between mechanical 
and optical degrees of freedom [I-Q. In particular, since 
the mechanical systems such as a dielectric membrane or 
sphere have masses much greater than that of an atom, 
the study of quantum optomechanics may test the foun- 
dation of quantum theory in macroscopic systems 161]. 
Typically, the systems considered in cavity quantum op- 
tomechanics are deformable cavities. These cavities are 
subject to radiation pressure pushing their cavity 'walls' 
apart, which in turn changes the field dynamics. For such 
systems, the deformation of the cavity is fundamental to 
the strong coupling between optics and mechanics. 

In this paper, we discuss another type of optomechan- 
ical coupling that results from the rotation of the optical 
cavity, in which case the coupling remains even in the ab- 
sence of cavity deformation. Physically, an optical field 
can affect the rotational motion via the electromagnetic 
torque exerting to a dielectric object [3-Q, and this has 
been studied in a sequence of experiments flOl4l4| . The 
mechanical rotation in turn affects light inside the dielec- 
tric, not only because the dielectric changes its orienta- 
tion, but also due to the motional-induced polarization 
and magnetization [ISj . Such a rotational optomechani- 
cal coupling could lead to a non-trivial coupled dynamics, 
and it is then a natural question on how the problem can 
be formulated self-consistently. In particular, a Hamilto- 
nian formalism of the system would allow a generalization 
to a fully quantized theory, in which both the optical and 
mechanical degrees of freedom are quantized. 

We remark that there are recent studies beginning to 
explore quantum effects in optically-trapped dielectric 
sphere [la - llSj . Since the orientation of the levitated di- 
electric particle is not fixed, rotational dynamics could be 
possible. However, it remains unclear about the strength 
of rotational coupling and what fundamental effects can 
be produced when light is acting on a rotating sphere 
and vice versa. 

The goal of this paper is to formulate a Hamiltonian 
that can address the rotational dynamics of the coupled 
dielectric-field system. Specifically, we consider a rigid 
dielectric sphere with radius R and moment of inertia /, 
placed in free space. The dielectric constant of the sphere 



is given by 



e(r) 



n 
1, 



|r-ro| <R 
otherwise. 



(1) 



We have used the convention eo = /io = 1 (i.e., c = 
1), and assumed non-magnetic dielectric /i = /iQ. We 
have also assumed a non-dispersive and non- absorptive 
dielectric. 

The sphere is free to rotate about any axis, but its 
center-of-mass (CM) is fixed at tq. In practice the sphere 
may be confined by an external potential so that the CM 
of the sphere moves about an equilibrium position. We 
assume such motion to be negligibly slow and of neg- 
ligible amplitude, then since for a spherical object, the 
CM motion does not directly couple with the rotational 
degrees of freedom [iSl], our approach here would be a 
good approximation. The system is specified by the La- 



grangian 
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where a; = (7 sin /3 cos a — (3 sin q;)x -|- (7 sin /3 sin a + 
/3cosa)y + {a + 7COS/3)z is the angular velocity of the 
sphere, a,/3,j are the 3 Euler angles specifying the 
orientation of the sphere (we follow the convention in 
Ref. [19,]). C is the Lagrangian density of the field 
after eliminating the electronic degrees of freedom of 
the dielectric. To find C, we go to an inertial frame 
5'(r) in which the dielectric element at r is instanta- 
neously at rest. Assuming the acceleration of the dielec- 
tric does not change its macroscopic properties, the field 
Lagrangian density at r in 5"(r) is given by the famil- 
iar form: C ^ ^ (eE'^ - B'^), where E' and B' are the 
electric and magnetic fields in S"(r), respectively. As the 
Lagrangian density is Lorentz-invariant, C can be read- 
ily obtained from the Lorentz transformation of the fields 
from 5'(r) to the laboratory frame S. We confine our- 
selves to a non-relativistic motion of the sphere, so that 
the velocity v(r) = u) x (r — rg) of the dielectric element 
at any point r satisfies |v(r)| <^ c. Linearizing C up to 
first order on a;, the Lagrangian reads 
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where T = / (fir{e — l)(r - Tq) x (E x B), which takes 
a form similar to the field angular momentmn stored in 
the sphere. The Lagrangian ([3]) is a generalization to 
that of Barton et al. 20] and Salamone [2l!] (in the case 
/x = 1), which consider a one-dimensional configuration 
and focus on CM motion of a dielectric slab. Here we 
will take a; as a degree of freedom which interacts with 
the field through the —us ■ T term. 

The electromagnetic field is specified by the scalar 
potential V{r,t) and vector potential A(r,i) under the 
generalized radiation gauge V • [e(r)A] = in the pres- 
ence of dielectric |22h24| . with E = — (dtA) — \7V and 
B = V X A. For the completeness of our theory, let us 
first discuss the Euler-Lagrange equations of the system 
before going to the Hamiltonian. 

The first Euler-Lagrange equation for the field is a re- 
statement of the Maxwell equation V • D = 0, which 
reads: V • (eVF) = V • [(e - 1) v x (V x A)]. This is 
understood from the fact that the polarization of a mov- 
ing dielectric element is P = (e — 1) (E + v x B) to first 
order oi v = |v|. Under the generalized radiation gauge, 

V is not a degree of freedom and it is determined by the 
instantaneous values of A(r,t) and v(r,t). We see that 

V is linear in w, and vanishes when the dielectric is at 
rest. The second Euler-Lagrange equation of the field is 
the wave equation: dt (edtA) -I- V x (V x A) = j, where 
j is a motion-induced (i.e. 0{v)) source current density. 



j = dt [-eVV + (e - 1) V X B] + V X M. 



(4) 



The wave equation is consistent with the Maxwell equa- 
tion V X B = dtT> + V X M in which M = -v x P is the 
magnetization of a moving dielectric with the polariza- 
tion P. From these two Euler Lagrange equations, we see 
that the motion-induced source terms are fundamental to 
the sphere-field coupling, without which the rotation of 
the dielectric sphere cannot affect the time evolution of 
the field. 

The mechanical equation of motion follows from the 
Euler-Lagrange equations of the Euler angles {dL/dC) = 
{d/dt){dL/d(), (^ = a, /3, 7. In terms of a;, we have 



rd^ ^ dT 

I-r = -u}xT+—-. 

dt dt 
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The two terms on the RHS of Eq. ([5]) characterize two 
different types of dynamics of ui. The first term de- 
scribes a precession about the T axis with a frequ ency 
|r|//, which keeps the magnitude of a; unchanged [25| . 
On the other hand, the second term may change the 
magnitude of U3 along the T axis. We remark that 
Eq. (O is consistent with the conservation of total an- 
gular momentum ^{Ii^ + 3 f) = which follows from 
the rotational invariance of the Lagrangian ([3]). To ze- 
roth order in u;, the total field angular momentum reads 
Jf = / rf^r(r - ro) x (E x B) [26f. 

We now turn to the Hamiltonian defined from L by 



H = 



E ^H + j 



n • {dtA) d^v - L, 
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where n(r, t) = [dC/d{dtA)] is the field canonical mo- 
mentum density, and p^ = (dL/dC) are the canonical 
momenta conjugate to the Euler angles. We introduce a 
canonical angular momentum J in terms of p^ [27]: 

J'x — ~ cot P cos apa — sin app + esc /? cos ap~f , 
Jy = — cot /? sin apa + cos apj^ + esc j3 sin ap^ , 
Jz = Pa- (7) 

Explicitly, n(r,<) and J are given by 

n = -eE- (e- l)(v X B) = -D, (8) 

J = lu-T. (9) 

Note that 11 is transverse as V-D = 0, and J differs from 
the kinetic angular momentum lu) for non-zero fields. 
The explicit expression of the Hamiltonian (jG)) reads 



^^ (j + r')' I f ,3 fn ^2 

H = - — —-- + - I d^ri + B^ 



(10) 



with r' given by 



r' = - / d^r 



e-1 



[(r-ro)x(nxB)]. (11) 



This Hamiltonian takes a form similar to the minimal- 
coupling Hamiltonian in electrodynamics, with T' some- 
how playing the role of vector potential in the kinetic 
energy term. We note that in writing Eq. ()10[) . we have 
neglected field-dependent terms that are quadratic in v. 
These terms resemble the kinetic energy of the sphere, 
and contribute to a correction of the moment of inertia 
/ due to the field. Such a correction is typically very 
small compared with / for fields well below the dielectric 
breakdown of the sphere. 

With the classical Hamiltonian (fTO|) . the canonical 
quantization of the system is readily achieved by pro- 
moting the dynamical variables C,, p^ (C = a,/3,7), A(r) 
and n(r) into operators by postulating the commutation 
relations: 



[A,(r),H,(r')]=zM?-(r,r') 



(12) 
(13) 



where (5|(r, r') is a generalized transverse S-function 
in the presence of dielectric [22|. From Eq. ([7]), the 
commutation relations of the Euler angles a, /?, 7 im- 
ply that J forms a quantum rigid rotor [l9|, which 
includes the angular momentum commutation relation 
[Ji,Jj] = ihcijkJk- The quantum Hanfiltonian takes 



the same expression as Eq. pI7|) . but with T' defined in 
Eq. pip symmetrized, i.e. with the bracketed term in the 
integrand replaced by [(r - Tq) x (11 x B — B x 11)] /2. 
In order to discuss field excitations in Fock 
space, we project the field operators onto a com- 
plete set of mode functions, namely the TE and 
TM mode functions (using spherical coordinates with 



C=a:/3,7 



origin at Tq) *;„(fc,r) 



.{E) 



{k,r)Xi, 



and 



*/„(fc,r) = (i/fc)V X [u'l'^\k,r)Xim{0,<l})] respectively, 
with X;m(0, (J)) being the vector spherical harmonics [2q . 

and M; {k,r) and wj {k,r) are radial functions sub- 
ject to appropriate boundary conditions across r = R: 
V • [e(r)*i™(ft,r)] = V • [e(r)*,^(fc, r)] = H^. Substi- 
tuting the normal-mode expansion into the Hamiltonian 
(we take h — \ from here on) , 



H = 



(j+n 

2/ 



Hf 



(14) 



where Hp = /rf^E/,™ Wfe[aL(fc)«i™(fc) + ^L(^)^;m(fc)] 
is the field Hamiltonian with constant terms removed, 
aim{k) and bim{k) are the annihilation operators for TE 
and TM mode photons with quantum numbers {k,l,m) 
and a frequency ujk = ck, respectively. 

Equations ((TU)) and ((T^ are main results of this paper. 
It is important to note that the form of F' [Eq. (ITTT) ] 
is very similar to the field angular momentum stored 
in the dielectric, apart from some proportionality con- 
stant. Therefore, approximately speaking, the first term 
of the Hamiltonian pUj) and p^ represents an angular 
momentum coupling, i.e., the interaction corresponds to 
an exchange of angular momenta between the field and 
the sphere. In particular, if the field is localized inside 
the sphere, we expect that T' should become a good ap- 
proximation to the total field angular momentum, up to 
a multiplicative constant. 

As an illustrative example of the rotational optome- 
chanical coupling, we apply the Hamiltonian ^1]) to 
a configuration in which photons occupy a whispering- 
gallery mode (WGM). In this case photons can be con- 
fined inside the dielectric cavity with a long life time 
due to multiple total internal reflections. For simplicity, 
we consider that the field excitation is dominantly con- 
tributed by TE mode photons with frequencies ujk ~ (^Oj 
where uiq = cko is a resonant frequency of a TE WGM. 
Assuming that such a WGM has a narrow line width Kc 
and LJo is well separated from all other TM mode fre- 
quencies, it is sufficient to include TE modes only in 
the Hamiltonian. Furthermore, since the optical qual- 
ity factor of the spherical cavity is typically very high, 
i.e. Q = A;o/kc ^ 1, it is instructive to consider the dy- 
namics of the system within a time scale short compared 
with K~^. In this regime, the leakage of WGM photons 
is negligible. Then the field Hamiltonian associated with 
TE WGMs with orbital quantum number I is given by 
lm=-! ^octiCm, where Cm is the cavity mode op- 



erator W 



dk 



ai, 
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k — kn 
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Here the index I for Cm is suppressed for compactness. 

With the help of cavity mode c™ operators, the F' 
operator contributed by the TE WGM is approximately 
given by. 



F'«A 
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= AS, 
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where L = — z(r — Fq) x V, Yir 
monies and 



are spherical har- 



A = TTKc 



R 



dr{e 



iy\u'i''\ko,r)f 



(17) 



is a dimensionless parameter determined by the mode 
amplitude inside the sphere. Numerical calculations with 
the parameters i? = 10 /im, n'^ = 2.31, I = 120, fco = 
27r/(743.25 nm) leads to A = 1.12. 

In writing Eq. ()16|) we have employed the rotating- 
wave approximation (RWA), so that fast oscillating terms 
such as cl^c]^, are dropped. However, these terms are 
responsible for photon generation in dynamical Casimir 
effect '31j|, and they should be retained if such a quan- 
tum effect become significant, for example, when u}(t) is 
rapidly changing with time. We also remark that the 
angular integral in Eq. (fT6| gives the selection rules for 
the rotational coupling. Together with [cm,c|^'] = 5mm', 
we see that S defined in Eq. (|16p satisfies the angular 
momentum commutation relations [S'i, Sj] — ieijkSk- 

By combining Eq. (fT4|) . the interaction with the TE 
WGM leads to a Hamiltonian (in a rotating frame where 
Hp is eliminated), 

i/. = A^i±^ + (l-A)g + A(A-l)|!. (18) 

which describes a coupling between the angular momenta 
J and S, with coupling strengths characterized by A. In 
addition, the Heisenberg equations of motion for S and 
uj (noting that J — luj — AS) become 

S = Aw X S and lev = A(A - 1) (w x S) (19) 

respectively, indicating that the optical (S) and mechan- 
ical angular momentum (/a;) precess about each other in 
the rotating frame. Moreover, for the coherent time scale 
concerned here (i.e., with negligible cavity field decay), 
the gallery mode photons cannot change the magnitude 
of o), and the sphere can only precess about the instan- 
taneous S-axis. 

An order-of-magnitude estimate of the mechanical pre- 
cession rate can be made by A(A — 1)(S)// w (n^ — 
l)Nhl/pR^, where N is the number of cavity photons, 
and p is the mass density of the sphere. Under the same 
numerical parameters above, and assuming N w 10^ be- 
fore dielectric breakdown, the precession rate would be 
on the order of 10~^ Hz. If the sphere spins coherently at 
a macroscopic rate and the WGM field inside the cavity 
is sufficiently weak such that /|(u;)| ^ |1 — A||(S)|, then 
Lj « (u)) behaves approximately as a constant (classical) 
vector. In this case the field dynamics described by S 
is governed by an effective Hamiltonian i?cff = ^{<^) ' S, 
which has the same form as that of a magnetic moment 
in an external magnetic field. Therefore photons initially 
occupying the cavity mode m experiences a Zeeman-type 
frequency shift of mK{ujz). We note that the frequency 
shift should be resolved from the cavity linewidth (i.e. 



A{u!z) ^ Kc) in order to be observable. With an optical 
field with quality factor of Q ~ 10^°, such a condition 
would require the dielectric sphere to spin at a rate of 
1 kHz. 

To conclude, we have established a non-relativistic 
Lagrangian and a Hamiltonian for a three-dimensional 
sphere-field system in the v <S^ c regime. By including 
the motion-induced polarization and magnetization pos- 
sessed by the dielectric sphere, we have self-consistently 
determined how a spinning dielectric sphere and quan- 
tized light fields are coupled. The sphere-field interaction 
is described by the Hamiltonian pUj) and (ITil) as a cou- 
pling of the canonical angular momentum J to the quan- 
tity r', which is proportional to the field angular mo- 
mentum stored in the dielectric sphere. We further illus- 
trated the rotational coupling using a WGM description 



of the Hamiltonian, and identified the constant A that de- 
termines the sphere-field coupling strength. Within the 
coherent time scale of the WGM photons, we have shown 
that the optical and mechanical angular momentum pre- 
cesses about each other, and the degenerate WGM mul- 
tiplets would experience a Zeeman-type splitting under 
a strong mechanical rotation. Although such effects are 
weak in general, they reveal fundamental features arising 
from the rotational degrees of freedom of the fields and 
the sphere. Our work here should provide a framework 
to further explore quantum phenomena and applications 
of such rotational optomechanical coupling. 
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